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The angular dependence of proton-proton correlation functions is studied in central ''°Ca-|-^''Ca 
and **Ca-|-**Ca nuclear reactions at E = 80 MeV/j4. Measurements were performed with the HiRA 
detector complemented by the 47r Array at NSCL. A striking angular dependence in the laboratory 
frame is found within p-p correlation functions for both systems that greatly exceeds the measured 
and expected isospin dependent difference between the neutron-rich and neutron-deficient systems. 
Sources measured at backward angles reflect the participant zone of the reaction, while much larger 
sources observed at forward angles reflect the expanding, fragmenting and evaporating projectile 
remnants. The decrease of the size of the source with increasing momentum is observed at backward 
angles while a weaker trend in the opposite direction is observed at forward angles. The results are 
compared to the theoretical calculations using the BUU transport model. 



I. INTRODUCTION 

The spectra of particles emitted in nuclear reactions 
can include contributions from a variety of dynamical 
and statistical mechanisms characterized by vastly differ- 
ent timescales. Dynamical emission typically occurs over 
timescales as short as 10~^^ seconds. Statistical emis- 
sion can extend to much longer times. The descriptions 
of dynamical and statistical emission mechanisms require 
completely different theoretical formalisms. This compli- 
cates theoretical interpretations of measured spectra, as 
most experimental observables do not allow a model inde- 
pendent distinction between earlier dynamical and later 
statistical emission. 

The correlation functions relevant to intensity interfer- 
ometry investigation [H-Qi however, do not suffer this 
limitation. Due to their ability to probe the space-time 
extent of the sources of emission, two particle correla- 
tion functions allow a distinction between early dynami- 
cal and later statistical emission. This has been used to 
probe the emission mechanisms of a variety of different 
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particle types for a number of different reactions studied 
over a wide range of collision energies [3-0] ■ 

The sensitivity of the two proton correlation functions 
to the space-time extent of the source arises from the 
mutual nuclear (attractive) and Coulomb (repulsive) in- 
teractions between the two protons and from the anti- 
symmetrized nature of their wave functions Q. Gates 
on the proton pair velocity provide information about 
the sources of these protons at different times during the 
reaction. Proton pairs with higher total momenta in the 
rest frame of the source preferentially reflect the space- 
time extent of that source at earlier emission times when 
the source is smaller. Smaller sources typically display 
larger and broader correlation functions [J|. In contrast, 
proton pairs with lower total momenta tend to be emitted 
at later times after the source has expanded and cooled. 
Such sources typically display narrower, and weaker cor- 
relation functions. Thus, correlation functions can track 
the time evolution of a cooling, expanding source. 

Transport models 0, 0] have revealed the existence of 
a sensitivity of two-nucleon (p-p , n-p and n-n ) correla- 
tion functions to the density dependence of the symme- 
try energy and some sensitivity to isospin in two-particle 
correlation functions have been observed [13. Physically, 
this sensitivity was shown to come from the effect of the 
symmetry energy on proton and neutron potentials and 
their influence on the emission times of particles dur- 
ing the pre-equilibrium stages of the collision [1]. This 
suggests that investigations of isospin effects on reaction 
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dynamics and their links to the density dependence of 
the symmetry energy 11|, may profit from a more clear 
understanding of the time characteristics of different par- 
ticle emission processes and by the capability of isolating 
emissions from the early pre-equilibrium stages of the re- 
action 

Stimulated by these ideas, we have measured p-p cor- 
relations over a wide angular and kinematic range with 
high statistics. In order to investigate the existence of 
isospin effects, we have compared results from reaction 
systems with different N/Z asymmetries, i.e. ''^°Ca-|-^°Ca 
(N/Z =1) and 48Ca+'^SCa (N/Z =1.4), at beam energies 
E/A=80 MeV. 

Studies have been performed for the ■^^Ar-|-'*^Sc reac- 
tion at E/A=80 MeV, providing some guidance for the 
dependence of the source on the momenta of the outgoing 
protons [l3j . The correlation functions measured in 
^^Ar-l-^'^^Sc reactions show a strong decrease in the source 
size with proton momentum, consistent with emission 
from an expanding and cooling participant source. BUU 
transport calculations generally reproduce these experi- 
mental trends. 

In our investigations of ^°Ca -|-^° Ca and *®Ca -I-*® Ca 
collisions, we have measured the correlation functions 
over a broader range of angles than previous measure- 
ments, and have studied in detail the momentum depen- 
dent two-proton correlation functions at different angles. 
The obtained results show that the applied momentum 
gates have strikingly different effects on the size of sources 
corresponding to particles emitted at forward angles as 
compared to those detected at backward angles in the 
laboratory frame. The measurements show a strong in- 
fluence of emission from the expanding, fragmenting and 
cooling spectator matter that was not evident in previous 
measurements. We also extract the fraction of protons 
emitted over short timescales during the collisions from 
the height of the correlation function and the integral of 
the imaged source distribution. We see surprisingly little 
sensitivity of these fractions to the angle or momentum 
of the measured protons. In order to better distinguish 
dynamical from statistical emission mechanisms, we also 
compared the extracted results to expectations of a BUU 
transport model [13]. 



II. TWO PARTICLE CORRELATIONS 

The two-proton correlation function probes the spa- 
cial and temporal information about the particle emitting 
source because the magnitude of the final-state interac- 
tions and anti-symmetrization effect depend on both the 
spacial separation and the relative momentum of parti- 
cles di. 

Theoretically, the correlation function is related to the 
space-time extent of the source by the angle-averaged 
Koonin-Pratt equation 0, [l3l 

C{q) = 1 + R{q) = 1 + 47r / A'(g, r)S{r)r'^dr, (1) 



where the two-particle source function, S(r), is the prob- 
ability of emitting two protons with a spatial separation 
r. In general protons are not emitted simultaneously. 
Then r in Eq. [T] refers to the separation at the time the 
second proton is emitted. The source function satisfies 
the following normalization condition 



47r / S{r)r'^dr = 1 



(2) 



The angle-averaged kernel K{q^ r) is given by 

i^(<?,r) = |</>,(r)|2-l. (3) 

where (f>q (r) is the two proton wave function measured at 
the separation distance between particles r and at pair 
relative momentum q, defined in the center of mass of 
the pair by: 



9 = M = \\Pi 



(4) 



Within this approach, the shape of the proton-proton 
correlation function is affected by the nature of the final- 
state interactions. The anti-correlation at low q is a re- 
sult of Coulomb repulsion. More importantly, there is a 
characteristic final ^Sg interaction peak at a relative mo- 
mentum of 20 MeV/c. If there are two distinct emission 
timescales, a fast dynamical emission and a slow statis- 
tical emission, the height of this final interaction peak 
should primarily reflect the fraction of protons emitted 
during the fast pre-equilibrium stage of the reaction and 
the size of the emitting source. The width of the peak at 
20 MeV/ c is solely affected by the space-time of the fast 
pre-equilibrium source fisj . Therefore a detailed study 
of the overall shape of the correlation function allows 
one to extract the space-time extent of the source and 
constrain the relative contributions from fast and slow 
proton emitting components. 

Experimentally the correlation function can be written 

as 



C(pi,p2)=AA 



B{pi,P2)' 



(5) 



Here, the numerator from Eq.[S]is the distribution of two 
protons with momentum pi and p2 detected in the same 
event. The denominator describes the uncorrelated back- 
ground distribution and is constructed using so-called 
event-mixing method (lol . [20| where each particle within 
a pair comes from a different event, taking into account 
the experimental two proton detection efficiency. J\f is 
a normalization factor, which typically results in corre- 
lation functions that are close to unity at large relative 
momenta HKH. 

We used two different methods to extract the sizes of 
the sources presented in this paper. In one method, we 
employed the imaging technique [23l - [25j to extract both 
the size of the source and the source distribution pro- 
file S{r) from the measured correlation functions. In the 
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other method, we obtained source sizes by fitting exper- 
imental correlation functions with the Koonin-Pratt for- 
mula (Eq. [1]) assuming the Gaussian source distribution 
S{r) given by 



Sir) 



X 



G 



(6) 



For the Gaussian source, there are three free parameters: 
the normalization of the correlation function, M, Ag and 
the source size Rq parameters of Eq. [6l 

According to Eq. [21 Ag ~ 1 if the emission of all 
protons used to construct the correlation function is de- 
scribed by the source function. While some protons are 
emitted over a short time scale after the collision and are 
strongly correlated, other protons can be emitted over 
very long timescales due to evaporation processes and 
secondary decays. Since the strength of the correlations 
reflects the spatial separation between the two protons at 
the time the second proton is emitted, early protons are 
not correlated with protons emitted at later times and 
late protons are only weakly correlated with each other. 

When both early and later emission occurs, the 
width of the peak in the p-p correlation function at 
20 MeV/cprimarily reflects the early emitted particles 
(fast source with smaller source sizes). Slowly emitted 
particles, coming from long-lived and more extended sec- 
ondary decay sources primarily influence the correlation 
function at low values [isj . If one is not primarily con- 
cerned with low q- values, these late emissions of protons 
largely reduce the magnitude of the correlation function 
while not usually strongly modifying its shape [l^- In 
this case. Equation [2] has a more general form that re- 
flects the fact that not all protons are correlated with 
each other, given to a good approximation by 



47r / S{r)r^dr = A. 



(7) 



The A parameter represents the fraction of pairs where 
both protons are emitted by the fast source represented 
by S{r) over the range of r represented in the integral in 
Eq. [7l The remainder 1 — A contains the contributions 
from pairs at large separation r outside of this range, 
where either one or both protons are emitted by the slow 
source at the late secondary decay stage of the reaction. 
The relevant proton pair fraction comes from the fast 
source; thus. A, can be well approximated by [Tsj 
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where / and the remainder 1 — / are the fractions of the 
total protons yields produced by the fast source and the 
slow source, respectively. 

To minimize apriori assumptions about the source 
function we follow the imaging techniques described in 
Ref. [3 I23l - l26| . and describe the source function S{r) 
by an expression involving three positive definite spline 
functions, which decreases monotonically with radius. 



We take the half-width half-maximum of the extracted 
source profile, ri/2j as a measure of the spatial extent of 
the source. This provides a simple size parameter that 
can be easily calculated even in the case of non-Gaussian 
source profiles where a Rq parameter as in Eq. |6] cannot 
be defined. In the specific case of a Gaussian source, the 
relation between r j /2 a-nd the size of the Gaussian source 
distribution (Eq. ^ is given by 



ri/2 = 2Vln2RG. 



(9) 



III. EXPERIMENTAL DETAILS 

We performed an experiment at the National Super- 
conducting Cyclotron Laboratory (NSCL), where pri- 
mary beams of '*"Ca (^^Ca) with an E/A=80 MeV im- 
pinged on ''^^Ca C^Ca) targets in a form of a thin mono- 
isotopic metallic foils of isotopic purity of about 97% 
(92.4%) by mass. We mounted the target near the center 
of the 47r Array, which housed 215 fast /slow phoswiches 
covering 85% of 47r solid angle around the target in the 
laboratory reference frame. The Att Array, with an in- 
side diameter of nearly 2 m, was instrumental in select- 
ing central events by requiring a high transverse energy, 
Et = Y.^E^sin^{^^) > 150 MeV [23. Here, 0, and E, 
correspond to the angles and energies of charged par- 
ticles detected in the 47r detector array. Assuming the 
transverse energy to monotonically decrease with impact 
parameter, this gate on Et corresponds approximately to 
an impact parameter range of < 6 {fm) < 4. 

In order to achieve both precise angular and energy 
measurements of the protons, required for correlation 
functions, we replaced one of the forward hexagonal mod- 
ules of the 47r Array with the High Resolution Array 
(HiRA) [ill. In our setup, HiRA consisted of 17 in- 
dividual telescopes in a hexagonal configuration, each 
housing a 65/iTO thin single-sided silicon strip detector 
followed by a 1.5 mm thick double-sided Si strip detec- 
tor with each face having 32 strips with a pitch of 2 mm 
and an active area of 62.3 x 62.3 mm^. The orthogonal 
orientation of the front and back strips of the thick Si 
detector, which was located 63 cm from the target, al- 
lowed for angular resolution of 59 « 0.05°. In order to 
allow the high-precision angular determination of the po- 
sition, we measured the position of the target and silicon 
strips in HiRA with the Laser Based Alignment System 
(LB AS) [2^. Additionally, the Si detectors were backed 
by a cluster of four 39 mm long CsI(TI) crystals which 
served as the calorimeters. For this paper, we only an- 
alyzed protons which stopped in the CsI(TI) crystals. 
This resulted in a proton momentum range of approxi- 
mately 200-500 MeV/c. The angular coverage of HiRA 
with respect to the beam was 18 < 9Lab (deg) < 58 in 
the laboratory frame and 30 < 0cm (deg) < 110 in the 
ccnter-of-mass frame. 




FIG. 1: Experimental correlation functions from 
(left) and **Ca+''*Ca (right). The upper panels include pro- 
tons with low total momentum of the pair (500-640 MeV/c) 
while the lower panels represent proton pairs with a high total 
momentum (740-900 MeV/c). The dashed-dotted lines repre- 
sent the results of the fit assuming the Gaussian source distri- 
bution. The solid lines are reconstructed correlation functions 
from imaging. The dashed lines represent the calculations as- 
suming the Gaussian source distribution with non-zero life- 
time (Eg. llOf) : see Sec. llVl for more details. 



IV. EXPERIMENTAL RESULTS 

The correlation functions measured in our experiment 
arc shown in Fig. [T] The left panels present results 
from '*°Ca+'*°Ca and the right panels from ^^Ca-F'^^Ca 
collisions. The upper and lower panels are for pro- 
tons with total momentum of the pair in the laboratory 
frame of 500-640 MeV/c and 740-900 MeV/c, respectively. 
The correlation functions at the most backward angles 
(33 — 58°) in the laboratory frame are represented by 
squares and at forward angles (18 — 26°) are shown as 
circles. The results at intermediate angles (26 — 33°) are 
plotted as diamonds. 

In order to get quantitative information about the pro- 
ton emitting source we use the imaging technique to ex- 
tract the imaged source function. The fits to the correla- 
tion function are shown as the solid lines in Fig. [TJ The 
corresponding extracted source functions are presented 
as the light cross-hatched and dark solid bands in Fig. [2] 
for 33 — 58° and 26 — 33°, respectively. In general, the cor- 
relation functions at backward angles have source func- 
tions that are larger and more localized around r = fm. 

Imaging allowed us to reconstruct source distributions 
only at backward and intermediate angles. The imag- 
ing technique fails at forward angles when the peak at 
q = 20 MeV/c is not well defined. If the source were 



FIG. 2: Comparison of imaging technique to Gaussian fit of 
p-p correlation functions for ''"Ca-f*'^Ca (left) and ***Ca+**Ca 
(right). The upper panels include proton pairs with low total 
momentum (500-640 MeV/ c) while the lower panels represent 
proton pairs with a high total momentum (740-900 MeV/c). 



Gaussian, the peak would become negligible for large val- 
ues of the Ra parameter in Eq. [HI e.g. Rg > 5 — 6 fm. 
Both the presence of sources with such large spatial ex- 
tensions, and large statistical errors in the correlation 
function make convergence of the imaging method diffi- 
cult to achieve at forward angles. 

The fit quality and the normalization of the recon- 
structed source distribution, Xj from Eq. [3 are given 
in Table [D If no constraints are placed on the shape 
of the source function, the imaging method can provide 
other solutions i.e. source functions S{r), with compara- 
bly small values of /dof , where dof w 30 is the number 
of data points minus the number of fit parameters. How- 
ever, some of those solutions have unphysical properties, 
such as S{r) < so we exclude them from the analysis 
and the error estimation of the A/ parameter. 

We also performed fits to the experimental correlation 
functions using Eq. [1] and assuming the Gaussian source 
function given by Eq. [S] The corresponding fits to the 
correlation functions are denoted by the dashed-dotted 
curves in Fig. [T] For the angular ranges of 6* = 26 — 33° 
and 6* = 33 — 58°, these fits are nearly indistinguishable 
from the fits obtained via the imaging procedure, the lat- 
ter shown as thick lines in Fig. [TJ In these fits there are 
three fitting parameters: 1) the size of the source, Rg] 
2) the \g parameter (from Eq. [H]); and 3) the normal- 
ization of the correlation function, M (from Eq. [5]) . The 
best fit parameters are presented in Table [J The source 
distributions obtained from the Gaussian fit are plotted 
as the solid, dashed-dotted and dashed lines for 18 — 26°, 
26 - 33° and 33 - 58° in Fig. H respectively. 
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TABLE I: Comparison of system size, angular and momentum dependence of results obtained from reconstructed source 
distribution with imaging method, Gaussian fitting procedure and BUU transport model simulations. 



The correlation functions reconstructed from imaging 
and obtained from the Gaussian fit are very similar and 
match the data well at most angles, as it is shown in 
Fig.[TJ For the lowest momentum gate at = 33—58°, the 
peaks in correlation functions for the Gaussian sources 
are narrower and their tails lie consistently below the 
data and the imaging results for g w 40 MeV/c. This 
gives rise to the slightly wider widths of the correspond- 
ing Gaussian sources shown in Fig. [5] for these data. For 
the other gates, the results for Gaussian and imaging 
analyses are very similar; in some cases, the source func- 
tions provided by imaging method are slightly more lo- 
calized at r = fm than the corresponding fits with 
the Gaussian source functions. At the most forward an- 
gles where the size of the source is large and the cor- 
relation effect is not as strong as in the experimental 
data collected at backward angles, it was not possible to 
constrain the source function adequately via the imaging 
technique. There we used the more constrained Gaussian 
source function in order to extract information about the 
space-time extent of the source. Fortunately, the similar- 
ity between Gaussian and imaging analyzes at the other 
angles provides support for us to use the Gaussian ap- 
proach and lends confidence to the information it pro- 
vides. 

To provide the simplified measure of the source, we 
characterize the extracted sources using (also used 
in e.g. [1^ m, H^) for each set of data and method used 
to extract the source distribution or its size. Results 
are presented in Table U for both reaction systems, both 
pair momentum ranges in the laboratory frame, and all 
three angular selections. With the exception of the lowest 



momentum gate at = 33 — 58°, the values for ri/2 are 
consistent between imaging and a Gaussian fit. 

The sources from the collisions with larger initial ge- 
ometry, "'^Ca+'^^Ca (N/Z=1.4), are systematically larger 
than those from '*"Ca±'^°Ca collisions (N/Z=l). The av- 
erage increase in source size with A somewhat exceeds 
A^/^, which suggests that the average freezeout density 
is somewhat lower for the ^^Ca than for "'*'Ca. Due to 
the large value of the neutron-proton cross section which 
significantly exceeds the p-p cross section, the relevant 
density for proton freezeout may be the neutron den- 
sity rather than the total nuclear density. In this case, 
the additional neutrons in the ^^Ca+^'^Ca may shift the 
freezeout to lower overall density. Such a shift refiects 
detailed differences in the transport of neutrons and pro- 
tons that could be used to extract information about the 
relevant in-medium cross section. Calculations indicate, 
however, that such effects are subtle and dwarfed by the 
qualitative difference between the sizes at forward and 
backward angles and we, therefore, defer such detailed 
model investigations to a latter publication. 

Clearly, the strong angular and momentum depen- 
dence of the extracted source size is a much more dra- 
matic trend. The observed large increase in the source 
size occurs at forward angles at velocities comparable to 
that of the beam. (The total momentum of two beam 
velocity protons is approximately 800 MeV/c.) Cor- 
relation functions of similar magnitudes have been re- 
ported for protons evaporated from heavy residues pro- 
duced in ^^^Xe+^^Al reactions [3l|, and for protons emit- 
ted at energies comparable to the Coulomb barrier in 
40Ar+i97Au [13 and i^SXe+^^^Au reactions [11]. In the 
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FIG. 3: ri/2 as a function of total momentum for *''Ca+''°Ca 
and **Ca+''*Ca collisions and all three angular ranges. The 
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by red closed circles while those from the Gaussian technique 
are shown as blue closed triangles. Source sizes from BUU 
are shown as black open circles. 



latter case, however, a relatively small fraction of fast 
protons (f=0.30) was reported. From the Xq and A/ 
parameters obtained from the Gaussian fit and imaging 
method we calculated the fraction / of short time scale 
emitted protons [l^ , according to Eq. |8l The results are 
summarized in Table HI We show that the values of this 
/ parameter are consistent between both methods. All 
proton fractions exceed 0.5 and there is very little mo- 
mentum or angular dependence of the / parameter. It is 
interesting that the large sources at 6* = 18—26° and total 
momenta of 740 — 900 MeV/c observed for both reactions 
also have relatively large fast fractions eg. / > 0.7. This 
implies that more than 70% of the two-proton emission 
occurs at relative separations of r < 15 fm. This appears 
to exclude significant contributions (>30%) from evap- 
orative emission at emission time delays much greater 
than about 150 fm/c. 

On shorter timcscalcs, the relative importance of a spa- 
tial expansion of the projectile remnant versus the effects 
of an extended lifetime cannot be distinguished with an- 
gular averaged correlation functions such as those pre- 
sented here. However, some information can be gleaned 
by considering the limits of an expansion followed by an 
instantaneous emission versus an emission that extends 



over timescale of the order lO's of fm/c. For the limit of 
instantaneous emission, we approximate the correspond- 
ing density by assuming that it is uniform with same 
RMS radius as the Gaussian distribution given by the 
best fit. In this approximation, the freezeout density 
would be approximately Qfre ~ Aspec/l^iV^Rc)^}- 



Assuming, the projectile contains As 
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prior to fragmentation and Rc = 3—9 fm, one obtains es- 
timates for Qfre of Qfre = 0.004— O.lpo- This is somewhat 
below the density range assumed by statistical simultane- 
ous multi- fragmentation models [3J, [3^ . It is also below 
the density range, p = 0.2 — 0.4poj extracted from d — a 
correlations for the participant source in ^^^Xe-|-^^''Au 
collisions [l^l • Both comparisons suggest that the source 
fragments over a non-zero timescale. 

Alternatively, wc assume that decay occurs from a 
spherical source with Rg and vary the timescale of the 
decay. Following Koonin [l^ we assume a Gaussian emis- 
sion time distribution: i.e. emission rate oc ea;p(— i^/r^). 
This leads to a source function of the form: 



Sir) 



Ar 



(4^)3/2i?2^^i?^ + 0.5(i;r)2 

Here, V = j^j, where 1^ —v —Vq is the magnitude of 
the velocity v of the center of mass of the two protons 
relative to the velocity Vq of the source, r± (r||) is the 
component of r perpendicular (parallel) to v. 

The beam momentum per nucleon is roughly equal 
to the average proton momentum for the data with 
18° < 6* < 26° and 740 < P (MeV/c) < 900. Thus, 
most of the protons within this gate must be preferen- 
tially emitted perpendicular to the beam leading to an 
estimated velocity of « 0.16c. In this scenario, the 
space-time extent measured for those particles is a com- 
bination of the spatial dimension (Rg) and the lifetime 
of the source (r) . The dashed lines in the lower panels of 
Fig. [T] correspond to the correlation functions obtained 
with the source distribution from Eq. I10[ where Rg = 3 
fm and r = 100 fm/c (w 3.3 x IQ-^^s) for -^OCa-h^^Ca 
and Rg = 3.5 fm and t = 135 fm/c 4.5 x IQ-^'^s) 
^^Ca-|-^®Ca reactions. 

Those calculations show reasonable agreement with the 
experimental correlation functions. This " lifetime" is rel- 
atively short for a statistical evaporation process, but 
comparable to the times for expansion and disassembly 
during a multifragmentation process [36l - f4l| . 

To illustrate the inconsistency of the large source 
sizes at forward angles with a straightforward dynami- 
cal origin, wc simulated '*°Ca-|-'*°Ca and ^^Ca-|-'*^Ca col- 
lisions at E/A=80 MeV. We chose a paramctrization of 
BUU such that an energy dependent in-medium nucleon- 
nucleon cross section reduction was employed [l^]. We 
also included momentum dependence in the mean field 
with a soft equation of state [1^. We chose the density 
dependence of the symmetry energy to be 7 = 0.7, which 
is in agreement with Rcf. 
duction of j4 <3 clusters 



441. We also included the pro- 
this tends to increase 
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''1/2 by approximately 1 fm. From the information pro- 
vided by the transport model we constructed the source 
functions for the same momentum and particle emission 
angle in the laboratory as used in the experimental anal- 
ysis. We included only the protons emitted at energies 
and angles that could have been detected in the exper- 
iment. We calculated the quantity ri/2 from the source 
distribution. We show the comparison between theoreti- 
cal and experimental values for ri/2 in Fig. [3l 

In Fig. [21 we can see that BUU can reproduce the ex- 
perimental data well at backward and intermediate an- 
gles for both pair momentum ranges measured in the 
laboratory frame, but underpredicts the sizes at forward 
angles, especially for protons in the high momentum gate. 
These high momentum particles move at close to the 
beam velocity. We have calculated source radii for a wide 
variety of different mean fields and nuclcon-nucleon cross 
sections, but have not been able to find a choice of trans- 
port parameters that result in significantly larger source 
radii at forward angles and beam velocity. Such large 
radii indicate emission from a source that is much larger 
or longer-lived or both compared to the source that can 
be predicted by a dynamical model such as the BUU 
approach. A long lived source could explain the discrep- 
ancy with the BUU calculations, however a very long- 
lived source is inconsistent with the large fast fractions 
/ > 0.7 deduced from our measurements. However, BUU 
suppresses many fluctuations that lead to rapid multi- 
fragment disassembles. The failure of the BUU to de- 
scribe the p-p correlations at forward angles and beam 
velocities provides a clear demonstration of the impor- 
tance of such processes in this kinematic domain. 

V. CONCLUSIONS 

We studied the angular and momentum dependence of 
p-p correlations for central '*"Ca+"'"Ca and "'^Ca-l-^^Ca 
nuclear reactions at E/A=80 MeV. We found a strong 
angular dependence within p-p correlation functions re- 
flecting the different space-time extent of the source se- 
lected. Sources observed at backward angles, in the labo- 
ratory frame, reflect the participant zone of the reaction, 
while much larger sources are seen at forward angles are 
dominated by expanding, fragmenting and evaporating 
projectile-like residues. The obtained results show a de- 
creasing source size with increasing momentum of the 
proton pair emitted at backward and intermediate an- 
gles. In contrast we observe a weak trend in the op- 
posite direction at forward angles. At some level, these 
trends are consistent. In the rest frames of the respec- 



tive sources, higher velocity protons arc more strongly 
correlated than their lower velocity counterparts, consis- 
tent with emission from expanding and cooling sources. 
The protons with small laboratory momenta at backward 
angles move slowly in the rest frame of the participant 
source. Protons with large laboratory momenta at back- 
ward angles move rapidly in the rest frame of the partici- 
pant source. In contrast, the highest momentum protons 
at forward laboratory angles are nearly at rest relative to 
the fragmenting projectile remnants, and the lower mo- 
mentum protons at forward angle are actually moving 
at a higher relative velocity to the fragmenting projectile 
remnants. In both angular domains, we therefore observe 
smaller sources for protons moving at higher velocities in 
the frame of the source. 

Long evaporation times are not consistent with the 
fast fractions extracted from the correlation functions 
forward angles. The time scales estimated from our 
correlation functions are consistent with bulk multi- 
fragmentation time scales that have been extracted by 
fragment-fragment correlation functions. 

Wc find that BUU transport calculations reproduce the 
data well at backward and intermediate angles, but un- 
dcrpredict to reproduce the source sizes at forward angles 
at high momentum. There the data are consistent with 
expansion, multi-fragmentation and subsequent evapora- 
tion. The failure of the BUU to reproduce the source 
functions for this case can be attributed to the suppres- 
sion of the fluctuations leading to multi-fragmentation in 
this approach. 

In all cases, the ^^Ca-(-'*®Ca reaction system results 
in larger sources than the "'"Ca-t-'^'Ca reaction system, 
which can be partly attributed to a sensitivity of the 
source distribution to the initial size of the projectile and 
target nuclei. However, the effect appears to be some- 
what larger than the A'^/'^ scaling expected from such 
geometrical arguments. 
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